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ABSTRACT 



In this paper we present two methods to evaluate non-factorizable corrections to pair-production 
of unstable particles. The methods are illustrated in detail for ly-pair-mediated four-fermion 
production. The results are valid a few widths above threshold, but not at threshold. One 
method uses the decomposition of n-point scalar functions for virtual and real photons, and 
can therefore be generalized to more complicated final states than four fermions. The other 
technique is an elaboration on a method known from the literature and serves as a useful check. 
Applications to other processes than IV-pair production are briefly mentioned. 
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1 Introduction 



With the start of LEP2, quantitative knowledge of the radiative corrections to the four-fermion 
production process e"'"e~ — > 4/ is needed The full calculation of all these corrections will 
be extremely involved and at present one relies on approximations [||], such as leading- log 
initial-state radiation and running couplings [^. Another approach is to exploit the fact that 
the corrections, in particular those associated with the production of an intermediate ly-boson 
pair, are important. This (charged- current) production mechanism dominates at LEP2 energies 
and determines the LEP2 sensitivity to the mass of the W boson and to the non-Abelian 
triple gauge-boson interactions. As such, one could approximate the complete set of radiative 
corrections by considering only the leading terms in an expansion around the W poles. The 
double-pole residues thus obtained could be viewed as a gauge-invariant definition of corrections 
to 'W-pair production". The sub-leading terms in this expansion are generically suppressed 
by powers of Tw/Mw, with Mw and Tw denoting the mass and width of the W boson. The 
quality of this double-pole approximation degrades in the vicinity of the VT-pair production 
threshold, but a few Tw above threshold it is already quite reliable 0. It is conceivable that in 
the near future a combination of the above-mentioned approximations will result in sufficiently 
accurate theoretical predictions for four-fermion production processes. 

In the double-pole approximation the complete set of first-order radiative corrections to 
the charged-current four-fermion processes can be divided into so-called factorizable and non- 
factorizable corrections 0], i.e. corrections that manifestly contain two resonant W propaga- 
tors as overall factors and those that do not. In view of gauge-invariance requirements, some 
care has to be taken with the precise definition of this split-up (see below). In the factorizable 
corrections one can distinguish between corrections to ly-pair production and W decay. In 
this paper we give a detailed account of the non-factorizable corrections that were used in the 
analysis of 0]. From the complete set of electroweak Feynman diagrams that contribute to 
the full 0{a) correction, we will therefore only consider the non-factorizable ones, both for 
virtual corrections and real-photon bremsstrahlung. To be more precise, since we are only in- 
terested in the double-pole terms we are led to consider only non-factorizable QED diagrams 
in the soft-photon limit. This means that we use simplified expressions for loop corrections 
and real-bremsstrahlung interferences, i.e. the photon momentum k'^ is neglected in the nu- 
merators and whenever possible fc^ is neglected in the denominators. It does not mean that 
the photon energy is taken to be small in the actual loop/phase-space integrations. In fact, 
for real bremsstrahlung the photon is treated inclusively and the energy is extended to infinity, 
which simplifies and approximates the phase-space integrals 0. The errors associated with 
this procedure are formally of higher order in the expansion in powers of Tw/Mw in the energy 
region where AE ^ M-^, since photons with an energy of 0{M\y) force the PF-bosons off the 
resonance Here AE stands for the distance in energy to the threshold. In the energy region 
where Vw <^ AE -C My/, the accuracy of the approximation becomes of 0{Tw I AE). Finally, 
near threshold, where AE k, Tw, our approximation breaks down. There the dominant correc- 
tion comes from the Coulomb effect, which was discussed in great detail in the literature (see 
e.g. [1,0,1). 
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Figure 1: Virtual diagrams contributing to the manifestly non-factorizable py-pair corrections 
in the purely leptonic case. The scalar functions corresponding to these diagrams are denoted 
by -Doi23! -D0124, and -^01234- 

The above approach is the same as the one adopted by the authors of Ref. who were 
the first to calculate non-factorizable VT-pair corrections. For the present calculations, we 
have used two different methods. One is an extension of the treatment in the other is 
a modification of the standard methods, which involves a combination of the decomposition 
of multipoint scalar functions and the Feynman-parameter technique. The results obtained 
with our two methods are in complete mutual agreement. However, in contrast to [Q, a clear 
separation between virtual and real photonic corrections has been made in both methods, which 
is essential to establish the cancellations of infrared and collinear divergences. This treatment 
reveals a significant difference between our results and those obtained by the authors of 0. 
Our final results do not contain any logarithmic terms involving the final-state fermion masses, 
whereas in the results of explicit logarithms of fermion-mass ratios occur (see discussion in 
Sect. 4.1 of Ref. P]). This difference can be traced back to the fact that although the fermion 
masses can formally be neglected in the absence of collinear divergences, they have to be 
introduced in intermediate results in order to regularize those divergences before dropping out 
from the final results. 

1.1 Gauge- invariant definition of non-factorizable corrections 

In order to give a gauge-invariant definition of the non-factorizable QED corrections in the 
soft limit, we first restrict ourselves to the manifestly non-factorizable corrections, i.e. those 
not having two resonant ly-propagators as explicit overall factors. In addition we will restrict 
ourselves to the simplest class of charged-current four-fermion processes, involving a purely 
leptonic final state: 

e+(gi)e-(g2) ^ + W-{p2) ^ u^{k[)t{h) + {k2)Mk'2)- (1) 

Whenever possible, all external fermions are taken to be massless. The relevant contributions 
consist of the final-final and intermediate-final state photonic interactions displayed in Fig. |l|. 
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In principle also manifestly non-factorizable vertex corrections exist, which arise when the 
photon in Fig. ^ does not originate from a PF-boson line but from the jWW/ZWW vertex 
(hidden in the central blob). Those contributions can be shown to vanish in the double- 
pole approximation, using power-counting arguments Also the manifestly non-factorizable 
initial-final state interference effects disappear in our approach. This happens upon adding 
virtual and real corrections, as will be briefly explained later. 

The double-pole contribution of the virtual corrections to the differential cross-section can 
be written in the form 



dcr^iit = 327rQ; Re 



i{p2 ■ ki)DiDoi23 + Kpi ■ k2)D2Doi2i + i{ki ■ A;2) ^1^2^01234 



dfJ-Born-, (2) 



where D12 = Pi2~ + iM^Vw are the inverse (Breit-Wigner) W^-boson propagators. The 
functions -D0123, -D0124, and -E01234 are the scalar integrals corresponding to the diagrams shown 
in Fig. |l|, with the integration measure defined as d'^k/{2'KY. The propagators occurring in 
these integrals are labelled according to: = photon, 1 = , 2 = W~ , 3 = and 4 = 
C.~ . Note that the factorization property exhibited in Eq. (§) is a direct consequence of the 
soft-photon approximation, which is inherent in our approach. As a result, the propagators 
hidden inside the central blobs of Fig. |l| are Born-like, i.e. unaffected by the presence of the 
non-factorizable photonic interactions. 

In a similar way, only interferences of the real-photon diagrams can give contributions to 
the manifestly non-factorizable corrections. The relevant interferences can be read off from 
Fig. |1] by taking the exchanged photon to be on-shell. The infrared divergences contained in 
the virtual corrections will cancel against those present in the corresponding bremsstrahlung 
interferences. 

The set of manifestly non-factorizable QED diagrams displayed in Fig. |l| is not U{1) gauge 
invariant, which can be explicitly seen from the non- vanishing of gauge-parameter-dependent 
terms when a general covariant gauge is used for the photon propagators. In order to achieve 
a gauge-invariant definition of the non-factorizable corrections, all (soft) photonic interactions 
between the positively (e+, ly"*", £"*") and negatively {e~ ,W~ charged particles should be 
taken into account, which also holds for hadronic final states. Looking at Fig. 0, this is equiv- 
alent to the set of all up-down QED interferences. The gauge invariance of these interference 
effects is a direct consequence of the fact that in the soft-photon approximation the matrix 
element can be viewed as a product of two conserved currents: one, where the soft photon is 
attached to the positively charged particles, and one, where it is attached to the negatively 
charged ones. Only three of those up-down QED interferences are already present in Fig. |l|, all 
others except one will vanish in the soft-photon approximation. In the soft-photon, double-pole 
approximation it is the "Coulomb" interaction between the off-shell W bosons that survives as 
an extra contribution to the differential cross-section (see Fig. ^: 



dcT^irt{Pi\P-2) = 327raRe 



i{pi ■p2)Ci 



012 



rfcTBorn- (3) 



The scalar three-point function Cou is defined according to the above-defined notation. The 
terminology "Coulomb" interaction should not lead to confusion. It is a contribution that is 
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Figure 2: The gauge-restoring "Coulomb" contribution. The corresponding scalar function is 



denoted by Coi2- In Sect.^]^ we shall briefly indicate the distinction between our "Coulomb" 
contribution, valid outside the threshold region, and the usual one, which is also valid inside 
that region. 



a part of the diagram in Fig. 0. In Sect. p.3| we shall briefly indicate the distinction between 
our "Coulomb" contribution, valid outside the threshold region, and the usual one, which is 
primarily valid inside that region. 

From the diagrams in Fig. |^ it is clear that we have to calculate four- and five-point scalar 
functions and related bremsstrahlung interference expressions, all in the soft-photon approxima- 
tion. In the next two sections we focus on the analytical results as obtained with the modified 
standard technique (MST), consisting of various elements. In particular the decomposition of 
the virtual and real five-point functions into a sum of four-point functions is given some special 
attention in Sect.|[[] As such, the basic building blocks of the MST are the four-point functions 
D and the related bremsstrahlung interference terms D^. A general relation between the two 
entities is discussed in Sect. ^]4[ As a final step we derive in Sect. |^ the relevant scalar four-point 
functions D in the soft-photon approximation by applying the Feynman-parameter technique. 
The related functions are obtained by using a "particle-pole" expression and performing 
certain substitutions. 

In Sect.^ we present the calculation along the lines of Ref. 0, which involves the method 
of direct momentum integration. This calculation will serve as a check of the MST results. 
Moreover, it is required for pinpointing the source of the differences with [Q. Whereas the 
method of Sects. ^ and |^ seems to be quite general, the method of Sect.| is unlikely to be 
applied to ra-point functions with > 5. As we will see, it just becomes too complicated. 

In Sect.^ we present the main analytical features of our study and we indicate how the 
results for semileptonic and hadronic final states can be obtained from the purely leptonic 
case. For the numerical implications of the non-factorizable corrections we refer to H]. Our 



"'^For higher n-point functions (n > 5) this decomposition can be carried out in an analogous way. Thus 



in principle the methods outlined in Sects. EJ^ and 2.2 provide the basic tools for considering more involved 



non-factorizable corrections, e.g. for six-fermion final states. 
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calculations confirm that non-factorizable corrections vanish in the special case of initial-final 
state interference, thereby making non-factorizable radiative corrections independent from the 
W production angle, and in all cases when the integrations over both invariant masses of 
the virtual W bosons are performed The practical consequence of the latter is that pure 
angular distributions are unaffected by non-factorizable 0{a) corrections. So, the studies of 
non-Abelian triple gauge-boson couplings at LEP2 [jlO| are not affected by these corrections. 
The non-factorizable 0{a) corrections, however, do affect the invariant-mass distributions {W 
line-shapes). These distributions play a crucial role in extracting the VT-boson mass from 
the data through direct reconstruction of the Breit-Wigner resonances. The non-factorizable 
corrections to the line-shapes turn out to be the same for quark and lepton final states, provided 
the integrations over the decay angles have been performed. Finally, in Sect.^ we draw some 
conclusions. 



2 Modified standard technique: basic ingredients 

In this section the basic ingredients are presented for the evaluation of non-factorizable correc- 
tions in the MST. As a first step we discuss the decomposition of virtual and real five-point 
functions into a sum of four-point functions. Subsequently we demonstrate how virtual and 
real contributions can be related in the soft-photon approximation. Having established the 
five-point decompositions and the relation between virtual and real contributions, the actual 
calculation in the MST boils down to the evaluation of scalar four-point integrals. 



2.1 Decomposition of the virtual five-point function 

In this subsection we derive the decomposition of the virtual scalar five-point function into a 
sum of scalar four-point functions. The derivation follows Ref. The reason for repeating 



this calculation lies in the fact that it will serve as guideline for the decomposition of the real 
five-point function, which has not been considered before. 

Let us consider the following general five-point function: 

77 _ f d^k 1 
i?oi234 - y NoN^N,NsN,' 

where 

iVo = A;^ - + io and Ni = {k + pif - m- + io. (5) 

Here io denotes an infinitesimal imaginary part. The plus sign accompanying this imaginary 
part is determined by causality. The mass parameter A is in principle arbitrary. In our case, 
however, it will denote a non-zero photon mass, needed for regularizing the infrared divergences. 

Before starting with the decomposition, we first derive a useful identity. To this end we 
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exploit Lorentz covariance and write 

The integral on the left-hand side is ultraviolet-finite and, when properly regularized, also 
infrared-finite. The quantities q on the right-hand side are therefore finite coefficients, depen- 
dent on masses and the invariants pi ■ pj = 1,2,3). Contracting this expression with the 
antisymmetric Levi-Civita tensor, one obtains the identity 



0, (7) 



(27r)4 N0N1N2NS 
where we introduced the widely-used notation 

In our convention, the Levi-Civita tensor is defined according to e'^^^s _ —5^-^23 = 1. The above 
identity will prove extremely useful in the following. 

The actual derivation of the decomposition formula starts with the Schouten identity 

aF = ^t;f (p.-fc), (9) 
1=1 

where the following notation was used: 

Ui — t , U2 — t , f 3 — fc , 6'4 — fc , u — cpjp2P3P4 ~ ■ y^^j 

Note that from the quantity a one can construct the Gram-determinant of the system, A4 = a^. 
The next step in the derivation is to contract the Schouten identity with k^, yielding 

4 

= J2{k-v,){k-pi). (11) 

i=l 

Now we can substitute 

A;2 = + - io, 
{k-p,) = ^[Ni-No-ri\, with r, =p2_^2^A2 (12) 

to arrive at ^ 

2a {No + X') = ■ -No- r,). (13) 

i=l 

In order to make the link to the scalar five-point function, one should divide this expression by 
N0N1N2NSN4 and perform the integration over d'^k. As a result of Eq. (|^ the Ni terms vanish. 
The terms J2{k ■ Vi) No can be transformed according to 



Y^[k ■ Vi) = e{fc+Pl){P2-Pl)(P3-pi){p4-pi) _ ,^14) 
i=l 



which can be verified by a direct check. After integration the first term will vanish. This can 
be most easily seen by making a change of integration variable, = +Pi, and subsequently 
applying Eq. (^. The complete expression now reads 

d^k aNo + 2aX'' + Y.ri{k-Vil^^ ^^^^ 



(27r)4 NQNiN2N^Ni 

The final step is to multiply this expression by a and to apply the Schouten identity and 
Eq. ([l^ to the last term in the numerator. This allows us to express the complete numerator 
in terms of the propagators appearing in the denominator: 



r d^k 2\''^^-\w^ + N^^^-\N^Y.{vi-w) + \Y.N,{vi-w) ^ 
J (27r)4 iVoiViA/'2A/'3A/'4 ' ^ ' 

with 



E^.<- (17) 
The final formula for the decomposition reads 



[w 



^ - 4A^A4) ^01234 = {W-Vi) D0234 + (W ■ V2) 1^0134 + (w ■ V^) L'oi24 + {w ■ V^) ^0123 



+ 



2A4 - ■ ^i) 



i=l 



D 



1234, 



where -D0234, -D0134, etc., denote four-point scalar functions containing the propagators with 
labels (0,2,3,4), (0,1,3,4), etc. 

The generalization of this decomposition to higher multipoint functions can be performed 
in a similar way [|ll|]. In general, a scalar iV-point function can be expressed in terms of the N 
underlying [N — l)-point functions. 



2.2 Decomposition of the real five-point function 

Using the derivation presented in the previous subsection as guideline, we can now try to derive 
a similar decomposition for the real five-point function. As can be read off from Fig. [1], by taking 
the exchanged photon to be on-shell, the real five-point function takes the form 

^^''''^ I i2nr2iu N[N^N!,N:,' ^^^^ 

where 



uj = \/k^ + X\ = Ni^2, and A^^^^ = Ar*^. (20) 

The photon is now on-shell, so k"^ = and Aq = 0. Note that the momenta pi, hidden inside 
Ni, are time-like and have positive energy components. 
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One can proceed in the same way as in the case of the decomposition of the virtual five- 
point function. The Schouten identity is still valid, but Eq. (|^ in its old form does not work 
in the case of real-photon radiation, and should be modified. In the derivation of the virtual 
decomposition, Eq. was used twice, leading to the nullification of 



(27r)4 N0N1N2NSN4 



and 



d^k No[E{k-Vi) + a] 



(27r)4 N0N1N2N3N4 



(21) 



In the case of real-photon radiation, this will correspond to 



d^k Y.{k ■ Vi) N[ 



{2-nf2uj N[N'2N'^Ni 



and 



d^k No [Eik ■ Vi 



(27r)3 2w N[N^NI^Ni 



(22) 



For the nullification of the second integral the validity of Eq. (^ is immaterial, since for the 
on-shell photon Nq = anyway. The first integral, however, is no longer necessarily zero. The 
fact that the photon is on-shell implies that k"^ = and that the propagators Ni are linear in k. 
By simple power counting, one can conclude that this integral is formally ultraviolet-divergent. 
For this reason, the Lorentz-covariance argument used in Eq. (^) is not correct any more and 
Eq. (|^) is invalidated. 

Apart from the modification of Eq. (|^) and the fact that A^o = 0, the derivation of the 
decomposition for the real five-point function is not changed, resulting in 



[w 



2 A\'A,)Ej' 



01234 



[w ■ vi) D^23A + ■ V2) /^S'134 + ■ V^) DI24 + (w ■ V^) D^, 



0123 



d^k Y.a{k-Vi)N[ 



(27r)3 2c<; N'-^NIiN'^N'^ 



(23) 



The main difference with the virtual decomposition is the occurrence of the last term in 
Eq. (PBI). It turns out that the poles in this particular integral can be moved in such a way that 
N!i Ni for all i. Indeed, the integral can be rewritten in the following way: 



d^k Y.{k-Vi)N[ 



(27r)3 2cj N[N'2N'^N'^ 



d^k 



(277)3 2W N1N2N3N4 



(24) 



with 



Ad) 
A(2) 



d^k {k ■ vi) Ni + {k ■ V2) N. 



(27r)3 2cu N1N2 

d^k ({k-V3), 1 
(27r)3 2u; j N1N2 




(25) 



Both A*^^^ and A*^^) are in fact zero. Let us consider, for example, one of the terms contributing 

(2) / d^k {k ■ Vi) 2iImA^3 



to A(2) = Af +Af ,e.g 



(27r)3 2cj A^iA'2 



N^N'^ 



(26) 
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This integral is ultraviolet-finite, even for an on-shell photon, and therefore no regularization 
is needed. Consequently, the Lorentz-covariance argument is valid: 



where the quantities q are finite coefficients. Contracting the last expression with v^^ = e^pipapsA" 
one arrives at A^^^'* = 0. Using similar arguments one can prove that A(i) = A(2) = 0. 

An important point in this line of reasoning was the use of Lorentz-invariance of the in- 
tegration d^k/uj. Such an integration is indeed Lorentz-invariant, provided that the area of 
integration is Lorentz-invariant. In the context of the double-pole approximation, the photon 
is treated inclusively, with the integration performed over all possible values of k up to infinity. 
If one would, however, consider an exclusive process, involving the introduction of a cutoff 
VLmaxi then the area of integration might fail to be Lorentz-invariant, and the decomposition 
stops at Eq. (|23|) . In order to successfully proceed beyond that point for exclusive processes, 
one should make sure that the cut-off prescription, which defines the area of integration, does 
not introduce new independent four-vectors in the integral. If this condition is satisfied, a 
generalization of the decomposition to exclusive bremsstrahlung processes should be feasible. 

So, in our approach, the following identity has been established: 

d^k Ejk ■ Vi) Nl ^ f d^k • Vi) 

{2TTy^2uj N[N^Nl^N^ J {27t)^2uj N1N2N3N4' ^ ^ 

As was already noted, the integral on the right-hand side is formally divergent. Its virtual 
analogue, being formally finite, vanishes: 

d'k nk■v^)N, ^ 

Performing a contour integration in the lower half of the complex A;o-plane, indicated by the 
integration contour Co in Fig. H, one can use this identity to relate the photon-pole contribution 
to the particle-pole contributions: 

^ . f d^k Y.{k-Vi)Ni f d^k _ , Eik-vAN^ , , 

-2m - — — ^ = -/- — - — Vole— (30) 

J {2tx)^2u NiN2N-iNi J {2-k)^Nq NiN2N^Ni ^ ' 

Here "Po/e" denotes the complex particle poles that should be taken into account. Note that 
the left-hand side of Eq. ( |30[ ) corresponds to the real-photon radiation integral that we are 
pursuing to evaluate. The integral on the right-hand side does not correspond to on-shell 
photons any more, since A^o 7^ 0. Now we can use the Schouten identity to substitute 

44 4 

Y.{k ■ Vi) Ni = J2ik ■ Vi) No + 2ae + ^ r, (A; ■ Vi) (31) 

i=l 2=1 i=l 

on the right-hand side. After some rearrangements one obtains 

. r d'^k Ejk ■ Vi) Nj _ r d^k j:{k-Vi)+2a r d^k ^ E ri {k ■ Vi) + 2aX^ 
{27r)^2uj N1N2N3N4 i (27r)4 A^iA'2A'3A'4 J{2ti)^No NiN2N3Ni 

(32) 
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photon particle 



Figure 3: Integration contour in the lower half of the complex kQ-plane leading to Eq. (|3 



The first integral on the right-hand side can be simplified with the help of Eqs. (|T^ and 
yielding 

2,./^i(^ = „f4!L_±__ + K, ,33) 



(27r)4 2cj NiN2N3N^ 



(27r)4 N1N2N3N4 



with 



7^ 



Vole ^ ^ 



In App.^ it is shown that the TZ term in (|33D, which consists of a combination of particle-pole 
contributions, vanishes. After that the decomposition for the real five-point function can be 
written in a compact form, analogous to the decomposition of the virtual five-point function 



01234 



0123 



+ 2iA4 



(27r)4 N1N2N3N4 



(35) 



Note that the last term on the right-hand side of Eq. ( P^D is exactly a virtual scalar four- 
point function with a coefficient 2iA4. In comparing Eq. ([T8| ) to Eq. (|35|) one observes certain 
similarities: the first four terms in Eq. (^) are the radiative analogues of their virtual coun- 
terparts in Eq. ([181) . One may naively think that the last term in (|35|) should not be there, 
since it does not correspond to photon radiation. In fact, it is a direct consequence of having 
ultraviolet-divergent integrals during the intermediate steps of the derivation. 

This concludes the derivation of the decomposition of the five-point function correspond- 
ing to inclusive bremsstrahlung. As was noted before, generalization to the case of exclusive 
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bremsstrahlung is possible, provided that the cut-off is introduced in such a way that no new 
independent four-vectors appear in the integrals. In analogy to what was remarked for the vir- 
tual decomposition, also the generalization to higher multipoint radiation functions is possible 
and rather straightforward. One should simply follow the approach of Ref. |Tl| for multipoint 
scalar functions. 



2.3 Application of the five-point decompositions 

We can now apply the five-point decompositions to the non-factorizable IV-pair corrections. 
The virtual scalar five-point function, corresponding to the third diagram in Fig. |1|, reads in 
the double-pole approximation 

-E01234 = 2A4 D1234 + {w -vi) D0234 + {w ■ V2) -D0134 + {w ■ V3) D0124 + {w ■ V4) -D0123, (36) 

with 

Vif, = +£pip2kif,, w'' = Div'^ + D2vl^, Ai = [ep^p^kik2?- (37) 

Comparison with Eq. ( plSf ) reveals that the terms —J^iw ■ Vi) -D1234 have been neglected, since 
they are formally of higher order in the expansion in powers of Tw/Mw- Note that the scalar 
four-point function -D1234 is purely a consequence of the decomposition (^). It does not involve 
the exchange of a photon and is therefore not affected by the soft-photon approximation. Since 
the factor 2A4/i(7^ is already doubly resonant, -D1234 should be calculated for on-shell W bosons. 
The scalar four-point functions -D0134 and i5o234 are infrared- divergent and should be calculated 
in the soft-photon approximation. 

In the same way, the five-point function corresponding to real-photon radiation is given by 

El,^, = 2tA, Df,,, + {w' ■ v[) D^,,, + iw' ■ v',) D^,, + {w' ■ iQ D^,,, + {w' ■ v',) D^,,,. (38) 

The four-vectors w' and v^ are defined as before, but for real-photon emission. This is equivalent 
to the following substitutions: pi —pi, ki —ki and D2 — > -Dg. The radiation function 
-^1234 is an artefact of the decomposition ( |38D and does not involve the exchange of a photon. 
It can be obtained from 1^1234 by the substitutions pi —pi and ki —ki. In the double-pole 
approximation, i.e. for on-shell W bosons, this implies the relation .0^34 = ilmDi234. This 
property ensures the cancellation of the virtual non-factorizable Di234-dependent corrections 
against the corresponding real-photon corrections, provided that the integration over the W- 
boson virtualities (i.e. -Di,2) is performed. This phenomenon is a general consequence of the 
soft-photon approximation p. 



2.4 Connection between virtual and real contributions 



At this point we have reduced the calculation of the non-factorizable corrections to the evalua- 
tion of virtual and real four-point functions. We can, however, go one step further and establish 
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a connection between the contribution from the photon-pole part of the virtual scalar functions, 
D"' , and the corresponding radiative interferences. To this end, we consider for example the 
contributions related to -D0123 and -D^i23- The contribution of the radiative interference to the 
cross-section (see Fig. ^ is given by 

w (D^ R /■ cPfc 327ra (pa • Di 

— * 

where fco = ^ = 1^1 • 

This has to be compared with the corresponding photon-pole part of the virtual correction. 
This contribution is evaluated in the lower half of the complex /cg-plane, where the photon pole 
is situated at kQ = u = \k\ — io: 

^ m7 o f dH: 327ra (pa ■ ki) Di 

rf^virt(/^oi23) - rf^BornKey j^-^ [D, - 2ip, ■ k)][D, + 2{p, ■ k)][-2{k, ■ k) + toY ^^"^ 

This can be rewritten in the form 

A (T^i \-A T) f 327ra {p2 ■ k{) DI 

da^irtiV,,,,) - daeorn Ke J - 2(j9i ■ k)^ + 2{p, ■ k)][-2{k, • k) - lo] ' ^^^^ 

Comparing Eqs. (|39|) and (^I]) , one can readily see that da-[ea\{D^i2z) can be obtained from the 
photon-pole contribution to the virtual correction (icTvirtl-Doias), by adding an overall minus 
sign and substituting]^ Di —DI. In a similar way c?(Treai(-D^]^24) can be obtained from the 
photon-pole contribution to the virtual correction (i(Tvirt(-Doi24), by adding an overall minus 
sign and substituting D2 — > —D2. The different substitution rule reflects the fact that we will 
determine -D0124 and -D0124 from -D0123 and -D0123 by substituting (pi, ki) ^ {p2, k2). Note that 
this is equivalent to evaluating -D0124 the upper half of the complex fco-plane. 

Also the "Coulomb" and five-point contributions can be treated in this way, bearing in mind 
that the coefficients of the five-point decomposition also depend on -Di,2- In conclusion, the 
following relation emerges. The radiative interferences can be obtained from Eqs. (0) and 
by adding a minus sign, by inserting the decomposition given in Eq. (|36D, and by substituting 



- in the -D0123, -D0134 terms: -D0123, -D0134 ^ -D0123, -Ddi34 followed by Di -D*^, 

- in the £'0124, ^0234 terms: /^oi24, ^0234 ^ ^0124, ^0234 followed by D2 -D^, 

- in the -D1234 terms: -D1234 — > -^^34 followed by D2 —D^, 

- in the C012 terms: C012 Cqi2 followed by Di —D\. 

Here both -D0124 and -D0234 are determined by substituting (pi, ki) ^ {p2, ^2) in the expressions 
for -D0123 and -D01345 respectively. As such, the above connection between real and virtual 
corrections implies that -Do^24 -^0234 evaluated in the upper half-plane. 

^Note that rfcTBorn is not affected by substitutions of the form Di —D* {i — 1,2), since it only depends 
on \D1D2?. 
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Note that the above-presented connection between the virtual and real non-factorizable 
corrections hinges on two things. First of all, the inclusive treatment of the bremsstrahlung 
photon, with the phase-space integration extending to infinity. Second, the fact that both 
virtual and real corrections are calculated in the soft-photon approximation, inherent in the 
double-pole approach. 

As mentioned before, manifestly non-factorizable initial-final state interference effects are 
also possible in our approach. As stated in Sect. |l.l| we will now briefly indicate why these 
effects vanish. Let us consider, for example, the initial-final state interference contribution 
corresponding to the photonic interaction between the positron [e"''(gi)] and the positively 
charged final-state lepton [^"'"(/ci)]. In the soft-photon approximation, the contribution of the 
virtual interference to the cross-section is 



dcrvirt(A./) = -(i(TBornRe J 



d^k 32i7Ta {qi ■ ki) Di 



(27r)4 [A;2 - + io] [Di - 2{pi ■ k)] [-2{ki ■ k) + io] [-2(gi ■ k) + io] ' 

(42) 

Note that all particle poles are situated in the upper half of the complex /co-plane. By closing 
the integration contour in the lower half-plane, one finds that the complete virtual correction 
is equal to the photon-pole contribution 

^ (D R f^Pfc 327ra {q^-k^)Di 

aa^irtl^./J a^Bornney ^^^^^^^ _ 2[p^ . k)][-2{k^ ■ k) + io\[-2{q^ ■ k) + to]' ^ ' 



with ko = u = y k"^ + — io. On the other hand, the corresponding bremsstrahlung interfer- 
ence can be written as 

, f d?k ?,2'na{qi- ki) DI 

da,,^,[D^,f) = -daBorn Re J ^^^^^^ [DJ + 2ip,.k)][2{k, ■ k) - io\\-2{q^ ■ k) + lo] ' ^ ^ 

By comparing the last two expressions, one can readily derive that the virtual and real interfer- 
ences only differ by an overall minus sign and the substitution Di — >■ —D\. In the next section 
we derive an explicit expression for infrared-divergent virtual scalar four-point functions [see 
Eg. ([7^)]. PI From this expression one can see that the substitution Di —D\ does not change 
the real part of the interference (^2]), i.e. the sum of virtual and real interferences gives rise 
to a vanishing non-factorizable correction. Analogously, no other non-factorizable initial-final 
and initial-intermediate state photonic interferences contribute to the double resonant cross- 
section, if both virtual and real corrections are included. Similar arguments can be used to 
prove that initial-state up-down QED interferences vanish in our approach. 



^In our example we need Eq. (72) with ^2 replaced by —qi. Note that, as a result of this substitution, the 
invariant Si2 becomes negative. 
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3 Modified standard technique: 
Feynman-parameter integrals 



In this section we present the calculation of the relevant virtual scalar functions, using Feynman- 
parameter integrals. In addition the photon-pole parts of these functions are given, from which 
the real-photon corrections can be extracted. The striking difference with the usual calculations 
of scalar integrals lies in the systematic application of the soft-photon approximation. 



3.1 Scalar four-point functions in the soft-photon approximation 

In this subsection we illustrate how to calculate a virtual scalar four-point function in the 
soft-photon approximation and how to extract the photon-pole part. Consider to this end 

J (27r)4 [A;2 + to] [2{p^ ■k) + D^ + lo] [2(p2 ■k)+D2 + lo] [2(^3 ■k) + D^ + ioy ^ ^ 

where Di = pf — Mf. In general, the energy components of the arbitrary momenta Pi are 
not necessarily positive. In contrast to the usual Feynman-parameter technique, where the 
Feynman-parameter transformation is applied to all propagators, we apply it only to propaga- 
tors that are linear in k: 

' rMO, (46) 







with 



(27r)4 [fc2 _ p + io][koE{0 - k- m + MO + 



10 



3 



The quantities A{^) and p^{C) are given by 

A(0=E6A and p^(0 = 2E6pr. (4J 



i=l i=l 



The energy component E{S,) of p^{C) can be positive or negative. However, there is a freedom 
to choose -E(0 ^ 0? because one can always perform a transformation of variables fco — > — ^o- 
Then 

r rl^h 1 

hUO = / 7^ ... ... , ' r (49) 



(27r)4 [fc2 _ p + io][-ko\E{0\ - k- m + MO + 



10 



3 



In the complex fco-plane the denominators give rise to poles. There are two photon poles, one 
in the upper and one in the lower half-plane. The second denominator gives rise to a "particle" 
pole in the upper half-plane, for any value of ^j. It should be noted that this combines the 
three particle poles present in (^5]), which could lie in the upper or lower half-plane. Closing 
the integration contour in the lower half-plane we get 

/virt(0 = / 7^ ^ , (50) 

J {2-nf 2\k\{-\k\\E{i)\-x\k\\p{i)\ + A{i) + ioY 
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where x = cos 6, with 6 being the angle between p(^) and k. It is not very difficult to perform 
the rest of the integrations in momentum space. The final result is 



^^'"^^^ A{o - ^om)\y ^^^^ 

As we have seen in the previous section, the real-photon radiative interferences can be 
obtained from the photon-pole parts of the virtual corrections. Let us therefore consider the 
photon-pole part of (|45|) in the lower half of the complex /cQ-plane: 



/}7 = -i I i (52) 

' i2n)^ 2\k\[2{p,-k) + D, + io][2{p2-k) + D2 + io][2{ps- k) + Ds + 10]' 



with ko = y k"^ — io. One can again proceed by introducing the Feynman parameters to obtain 

= / rSi ^7 h ^ ^ Vr ^53) 

^^''> 2\k\[E{i)^k^-io-x\k\\p{0\ + A{0+^o\ 

Equations and (^) are the same up to small modifications. In the case of the full virtual 
scalar function, Eq. ( ^0]) was obtained after contour integration in the complex fco-plane. In 
that case we had the freedom to choose the contour in such a way that ^ 0. Now we have 
no such freedom. So, E{^) cannot be considered as a negative quantity any more. It is clear 
that the final answer will be 

i 1 

n-Ao = Aiom)+^oE{or ^^^^ 

This expression is very similar to the one derived for the full virtual scalar function. It can in 
fact be rewritten as 

where the ffist term in the curly brackets corresponds to the full virtual scalar function and the 
second term is the necessary modification. The second term in Eq. (|55D is the analogue of the 
"particle" -pole contribution in the approach of 0. Note that this term has an extra factor i. If 
all quantities were to be real (stable-particle case), then this term would not contribute to the 
non-factorizable correction to the cross-section, for which only the real part is important. In 
the case of unstable particles, this "particle" -pole contribution is felt by the imaginary parts of 
the VT-boson propagators, resulting in a potentially non-zero contribution to the cross-section. 
If one were to evaluate the photon-pole part of (^51) in the upper half of the complex /co-plane, 
one merely would have to replace E{^) by —E{C,) in Eqs. ( ^4|) and (|55|). 



In practice, we calculate the relevant four-point functions -Dvirt as well as the corresponding 
particle-pole contributions -Dvirt*- The photon-pole part -D^irt is obtained as -0^;^ = -Dvirt— -0^1^*; 
which can then be used to evaluate the real-photon radiative interferences. The complex half- 
plane where the particle-pole (photon-pole) contributions should be evaluated is fixed according 
to the rules given in Sect.OI. 
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3.2 Calculation of the virtual scalar four-point functions 



In this subsection we present the calculation of the virtual scalar four-point functions and the 
associated photon-pole parts. Everything is considered in the soft-photon approximation. Since 
the four-point function -D1234 does not involve this approximation, we defer the corresponding 



results to App. [A.l| and merely refer to the literature 0, |l^ for its derivation. 



Before listing the various results, we define our notation. To write down the analytical 
results we need to introduce some kinematic invariants: 

"^1,2 = ^1,2) S = (]9l+P2)^ Si2 = {ki + k2f, 

S211' = {k2 + ki + k[)\ su2' = {ki + k2 + k'^f, (56) 
and some short-hand notations: 

yo = Y^, Xs = ^—^ + io, (3= a/i -4M^/s , 
U2 p + i 

/■ 1 '^122' . -1 ■5211' . 

3.2.1 The virtual infrared- finite four-point function 

We start off with the calculation of the infrared-finite scalar four-point function -D0123, which 
corresponds to the first diagram shown in Fig. |I[ This function is infrared-finite owing to the 
presence of finite decay widths in the propagators of the unstable W bosons. In the soft-photon 
limit we find 

~ J (2^ [A;2 + ^o] [Di - 2(pi ■ k)] [D2 + 2{p2 ■ k)] [-2{k^ ■ k) + lo] ' ^^^^ 

where -Di^2 = Pi 2 ~ M^^ + io. Originally the quantities D12 are real, with the usual infinitesimal 
imaginary part. At the end of the calculation the analytical continuation to finite imaginary 
parts can be performed. Then D12 = P12 ~ + iMw^w- 



Applying the Feynman-parameter technique as explained in Sect. |3.1| , we obtain the follow- 
ing representation 

with 

A(o = 6A+6^2, p'^(o = -2eirf + 26p^-2e3fcr. (6o) 

As was indicated before, the integral will be calculated for small final-state fermion masses 
and in the double-pole approximation. This implies 

2(pi ■ /ci) ~ ~ and 2(^2 ■ A;2) ^ P2 ~ (61) 
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$2(s+Si22')— Sl22' 



S S-i nnf 




Figure 4: The integration area (shaded region) in the (^,^2) Feynman-parameter space for the 
calculation of the particle-pole part -D0123 of the infrared-finite scalar four-point function i5oi23- 



The thicker curves indicate the solutions of p (^) = 0, with ^2 



/Sl22'- 



What is left is the integral over the space of Feynman parameters. The details of the 
integration are presented in App. |A.^ . The final answer reads 



0123 



2Ch 



£i2 x,;C +/:i2 -;C + 



1 1 

Z/o' C 

In 



Ml 

2 

mf 



1/0 

+ 21n(C) 



- >Ci2 — ; — 

\Xs yo 

ln(yo)+ln(C) 



The function Ci2{x;y) is the continued dilogarithm 

Ci2{x; y) = Li2(l - xy) + ln(l - xy) [\n{xy) - \n{x) - ln(y)], 



(62) 



(63) 



with Li2(a;) the usual dilogarithm and x,y lying on the first Riemann sheet. The answer for 
the second infrared- finite scalar four-point function, -D0124, can be obtained from Eq. (|62D by 
substituting {pi.h) ^ {p2,k2)- 



3.2.2 Photon-pole part of the infrared-finite four-point function 

As was explained in Sect.|3]l|, the photon-pole part of a scalar function can be obtained from 
the full scalar function by subtracting the "particle" -pole contributions. According to Eq. (0), 
the "particle" -pole contributions in the lower half of the complex ko-plane are given by 

D'ot^s = ^ / ^ Om)] S\p\a (64) 
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The integration area is defined by the ^-function for the energy and by the condition ^1 + ^2 < 1- 
The allowed area of integration and the curve where the 5-function has a non-zero value are 
schematically shown in Fig.^. The depicted situation represents the most general case for the 
kinematics we are interested in. 



After the integration over the (5-function has been performed, one is left with a simple 
one-dimensional integration of logarithmic type. The final result is 



7-) part 
-^^0123 



1 



8nM^ D2 - CDi 



ln(l - uoXs) - ln(l - Xs/Q 



(65) 



3.2.3 The virtual infrared-divergent four-point function 

In this subsection we describe the calculation of the infrared-divergent scalar four-point function 
-D0134, which enters the non-factorizable corrections through the decomposition of the five-point 
function. Similar four-point functions may also appear in the initial-final state interactions, 
but, as was mentioned before, the interactions of this type vanish in the sum of virtual and real 
contributions. The four-point function 

~ J (2^ [A;2 - A2 + 10] [Di - 2(pi • k)] [-2(A;i ■ k) + 10] [2{k2 ■ k) + 10] ' 

is infrared- divergent, because only one unstable particle is involved, which is not enough to 
regularize the divergence. Therefore we introduce a regulator mass A for the photon in order 
to trace the cancellation of infrared divergences in virtual and real corrections. Again we can 
apply the Feynman-parameter technique as explained in Sect. |3.1| . However, special care has to 
be taken with the photon mass A. As usual we can introduce Feynman parameters according 
to 

1 3 

D0134 = 2 j d'idl^l - ^e.) /oi34(0, (67) 


/o,34(^) = / 1 ^ (68) 

J {2'KY[k2-\2 + ,o\[-ko\E{i)\-k-p{i) + A{i)+iof 

and 

^(0 = 6^1, p'^(0 = -2ei^r-26rf + 2e3A;^ (69) 

Again we can exploit the freedom to perform the variable transformation /cq — > — /cq in order 
to fix the sign of the energy component E{^). After the integration over momentum space, the 
details of which can be found in App. |A.3| , we obtain 

d { 1 , /A- V^^-AV 



with 



with 



/oi34(0 - ^ i vA^-AV ^"" U + V^^-AV J i 

^ = ilMl, + iWl + ^3^2 + - 66^12 + 66(M^ - S21lO. (71) 
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Figure 5: The area of integration in the (M,t)-plane for the calculation of the particle pole 
-^om- The shaded region is the area of integration where E{u,t) > 0. The doubly-shaded 
region is the area of integration where E{u, t) > and p'^{u, t) > 0. The thicker curves indicate 
the solutions of p^{u,t) = 0, with u* = — ('M^/m^. 

The masses of the final-state fermions, mi and m2, are taken to be small in our approximation. 
One is left with a twofold Feynman-parameter integration (see App. [A.3|) , which results in 



0134 



1671^512 Di 



Li2 1 + 



Sl2 



21„l^lh> 

m2 



17111712 

-Si2 - io 



71^ , 2fMw\ , 2 

H h n^ — — + n^ 

3 V mi y KQ'Mw 



(72) 



The answer for the second infrared-divergent scalar four-point function, -D0234, can be obtained 
from Eq. ( [T^D by substituting (pi, ki) ^ (^2,^2)- 

3.2.4 Photon-pole part of the infrared-divergent four-point function 



The procedure for calculating the photon-pole part -D0134 follows our general strategy, i.e. we 
calculate the corresponding "particle" -pole contributions and subtract them from the full virtual 



scalar function. The difference between this case and the one discussed in Sect. 3.2.2 lies in the 



fact that here both photon-pole and "particle" -pole contributions are infrared-divergent. So 
as to keep track of the cancellation of the infrared divergences, we again introduce a regulator 
mass A for the photon. 
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In App. |A.3| the following convenient representation for /oi34(0 was derived in Eq. (|155|) : 

dz 



'0134 



(0 



-\E{0\ v/i^+A^ - I m \z + A{0+to 



2 ■ 



(73) 



From the general discussion in Sect. |3.1| we know that the photon-pole contribution to the four- 
point functions in the Feynman-parameter representation, I^^^^{C,), can be obtained from the 
complete virtual function, IvhtiO^ by a substitution ~^ ^-^(0- Then the "particle" -pole 



contribution, DqiI\ = -D0134 ~ -Dqi 



has the form 
1 



(74) 



where 



E{i) dz 



(75) 



By analysing the pole structure of the last expression, it is easy to realize that /oi34(0 is only 
non- vanishing for p^(^) > 0. This results in a simple formula 



em)] d 



4. ^.^(0 h(o,A 



[A(0+io]- 



(76) 



What is left is the integration over the space of Feynman parameters. In order to simplify the 
calculation it is advisable to make the change of variables ^1 = t/{l+t+u) and ,^2 = I / (l+t+u) . 
The area of integration in the {u, t)-plane is shown schematically in Fig. ^. The final result is 



7-) part 
-'^0134 



87rsi2-Di 



w 



In 



(-) 

\m2/ 



(77) 



This is the same result as obtained in Sect.H, where we will approach the various calculations 
in a significantly different way. 



3.3 The Coulomb-like scalar three-point function 

As was pointed out in Sect. |l.l| , a gauge-invariant definition of the non-factorizable corrections 
requires the proper inclusion of a Coulomb-like contribution. In this subsection we calculate 
the associated scalar three-point function in the soft-photon approximation. This three-point 
function is infrared- finite, but ultraviolet-divergent. This divergence occurs as a result of the 
fact that we neglect the dependence of the propagators, following our general soft-photon 
strategy. Although the virtual and real Coulomb-like contributions to the cross-section are 
separately ultraviolet-divergent, the sum is finite. 



20 



The virtual Coulomb-like scalar three-point function C012 is defined as (see Fig. 0): 



012 



(27r)4 [P + io\[D, - 2(pi • k)][D^ + 2{p^ ■ k)] ' 



(78) 



Similar to the calculation of the scalar four-point functions, we introduce the Feynman param- 
eters only for the propagators that are linear in k. We limit the area of integration over k by 
the condition \k\ < A, where A ^ F^y. After having performed the momentum integration, 
one is left with a one-dimensional integral over Feynman parameters: 



012 



1 



\p\{\E\ + \p\) 



In 



p\{\E 
E\ + \p 



\P\) 



In 



\E\ - \ p\ - A/ A 
-A/A 



A/A_\ 



(79) 



where A{^) = 6 ^1 + 6^2 and p^(0 
Feynman parameters yields 



-A/A 

—2C,iPi + 2^2 P2- The final integration over the 



012 



167r2s/5 



>Ci2U/o; 



Xq 



In 



-iDi 
2MwA 



{-•- 



-2Li2 1- 



Xq 



+ In 



-tD2 
2MwA 



— 2i7r In 



In'(yo) 



l + Xs 



10) 



In a similar way one can calculate the particle-pole contribution Cqi2^, which can be used 
to extract the photon-pole part C012 = C012 — CqS*- The final answer for the particle-pole 
contribution reads 

C^i2 = ^{1^(1 - ^s) + ln(l + X.) - ln(l - Vox,) - In(^) }• (81) 

Note that the real part of 6*0^2* does not contribute to the non-factorizable corrections. There- 
fore, the A dependence effectively drops out from the particle-pole contribution. 

One may wonder how the non-factorizable contributions (0) and (|^) to the cross-section 
compare with the Coulomb contribution as calculated in the literature. Our calculation is 
based on the assumption of being at least a few widths away from the threshold [the accuracy 
of this approximation is of 0(r^/A£')], whereas the Coulomb effect in the literature is valid 
in the non-relativistic region, where it approximates the cross-section with accuracy 0{l3). 
Therefore one could try to compare them in an overlapping region, where Tyy <C AE <C M^y. 
The Coulomb effect calculated in Ref. [§] consists of two parts. One contribution that is also 
present for on-shell VT-bosons, and one that comes from the off-shellness. The former is related 
to factorizable corrections and the latter is related to non-factorizable corrections and vanishes 
upon integration over the virtualitie^. In the overlap region {Tw AE -C Mw) this off-shell 
term equals the 1//3 part of our full expression. 

^Only soft virtual photons contribute to both parts. For the on-shell (factorizable) part of the Coulomb 
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4 Direct momentum-integration method 



In this section we present an alternative method of calculating the non-factorizable corrections, 
i.e. the so-called direct momentum-integration (DMI) method. As in Sects. |^ and ^ we use the 
soft-photon approximation and assume the charged final-state fermions to be massless, which 
is a good approximation for the process under consideration. In contrast to Sects. ^ and ^ we 
do not make any assumptions about the mass of the final-state neutrinos, because it does not 
simplify the calculation significantly. This gives us the opportunity to apply the results of this 
section to top-quark pair production. Following the approach of 0, we write the amplitudes 
corresponding to the diagrams shown in Fig. |I| in terms of virtual and real scalar four-/five- 
point functions. In contrast to the Feynman-parameter approach of Sect.^, we do not introduce 
Feynman parameters, but perform instead a direct integration over momentum space The 
calculation can be considerably simplified by an appropriate choice of the frame. 

First we calculate the infrared-finite virtual and real four-point functions. The calculation 
is close to the one presented in |Q, but in contrast to we make a clear separation between 
the virtual and real contributions. Our final result agrees with the one of as well as with 
the one obtained in the MST in Sects. 3.2.1 and p.2.2| . Next we calculate the infrared- divergent 



virtual and real four-point functions. Again we perform a separation of the real and virtual 
contributions, and provide a careful treatment of the divergences. All this is needed in order to 
trace the cancellations of infrared and collinear divergences. We find complete agreement with 
our results obtained in Sects. [3l273| and |3.2.4 At the same time the structure of the divergences 



in our results appears to be significantly different from the one obtained by using the method 
of [^, even in the complete answer when virtual and real corrections are summed up. Although 
the infrared-divergent scalar four-point functions do not appear directly in the answer for the 
non-factorizable corrections, the observed difference with the method of turns out to be 
indicative, because similar problems arise in the evaluation of the (infrared- divergent) five- 
point functions. Finally we calculate in the same way the virtual and real five-point functions. 
After summing up the virtual and real five-point functions, we find, in contrast to the result in 
0, that all collinear divergences cancel exactly, even for cross-sections that are exclusive with 
respect to the virtualities of the W bosons. 

In conclusion, the calculation presented in this section is an extension of the method of 
0. We provide a proper treatment of the infrared and collinear divergences, and make a clear 
separation of the virtual and real corrections. Because of this, the calculation becomes much 
more involved. We use the results obtained in this section as an independent check of the results 
of Sects. @ and |^. Although the methods are completely different and the answer of this section 
is very complicated, a perfect numerical agreement between our two calculations is observed. 

Before listing the various results, we first define the notation. For the calculations in the 



effect photons with momenta uj « P'^AIw and \k\ « /3Mw are important (hence can not be neglected 
in the propagators of the unstable particles). On the other hand, only photons with momenta to ~ Tw and 
|fc| « T\y/P give the leading contribution to the off-shell part of the Coulomb effect. Far from the threshold, 
where Tw ^ <C Mw, the two regions in the photon- momentum space are well separated. Because of this 
the effects are additive. Near threshold, where Tw ~ ^E, the two regions start to intersect. 
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DMI method we need to specify the momenta in the centre-of-mass frame of the initial state. 
Because of the soft-photon, double-pole (-Di,2 ^w) approximation, the four- momenta of the 
two intermediate W bosons are related in a simple way: 

p>t = {E,p) = E{l,v) and = (E, -p) = E (1, -^7), (82) 

with \v\ = [3 the (on-shell) velocity of the W bosons [see also Eq. (|57|)]. The other relevant 
momenta are 

k>t = {E^,h) = Ei{l,Vi) and = (^2, fca) = ^2 (1, ^^2), (83) 

with \vi\ = Vi = y^l — mf / Ef for {i = 1, 2). In addition we need the definition of some (polar) 
angles with respect to the direction of the boson: 9i = L{v, Vi) and Xi = cos 6*4 for {i = 1, 2). 
The difference of the azimuthal angles of ki and ^2 is given by 0i2. So, for sin 012 = the final- 
state three-momenta ki and k'^ lie in one plane. In the plane spanned by ki and k2 we define 
612 = l{vi,V2) and = cos6'i2. 



4.1 Non-factorizable infrared-finite corrections 

In this subsection we briefly describe the calculation of the infrared-finite four-point functions 
in the DMI scheme, following Ref. The result agrees with the one presented in 0, so this 
subsection is merely presented for completeness. The contribution of the infrared-finite virtual 
four-point function Z)oi24 to the non-factorizable matrix element is given by 

tTj ^ [ ^'^^ IGvra; (pi ■ /C2) (aA\ 

""Wj (2^ + to][2{k ■ k2) + io\[D^ - 2(pi ■ k)]{D2 + 2(p2 ■ k)] ' ^ ' 

where Mb = Mb/ {D1D2) is the Born matrix element of the process, involving the produc- 
tion of an intermediate 14^-boson pair and its subsequent decay. We start the calculation by 
decomposing the unstable W^-boson propagators according to 



1 



[D,-2{p^-k)][D2 + 2{p2-k)] 



1 

+ 



D^-2{p^-k) D2 + 2{p2-k) 



^ (85) 



D + ip-k 



where D = Di + D2- The first term has two particle poles: one in the lower and one in the 
upper half of the complex fco-plane. We close the contour in the lower half-plane, resulting in 
one particle-pole and one photon-pole contribution. The second term has all its particle poles 
in the lower half-plane. By closing the integration contour in the upper half-plane, only one 
of the photon poles will contribute. Note that the above decomposition mixes photon- and 
particle-pole contributions. In order to avoid possible confusion with the pure photon- and 
particle-pole contributions, we will write M0124 and M0124 if the decomposition is used. 
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4.1.1 Particle-pole residue 



We first concentrate on the particle-pole residue contributing to the first term in Eq. (pS]). This 
particle pole is situated at fco = V2-k. The corresponding residue reads 

r'part' 47ra 1 — (3x2 f 1 



Mo?24 = Mb ^ / ^ — ^- ^. (86) 

The propagators can be exponentiated by introducing an integration over "time" : 



Ana I - 13x2 7, , r S'k 



e 



ik-f 



Mjlll' = -Mb ^ / rfr rfn eH#-+2i-iJ / -1^ ; ^, (87) 

where 

f=2Tv~Ti{v2-v). (88) 

The integral is infrared- finite, so there is no need to introduce a non-zero photon mass as infrared 
regulator. The spatial integration can be recognized as a relativistic Coulomb potential of a 
moving particle: 

</,(r) = -47r / — — —y- ^ = = (89) 



Here ry and r±_ are the absolute values of the components of r parallel and perpendicular to V2'- 

r\\ = 2(3x2T - {I - (3x2)ti, = /3(2r + n) sin ^2. (90) 

Note that 1 — f 2 = m\/ E2 is small, but finite. 

To do the remaining integrations over r and ri, we can make a change of variables according 
to (t, Ti) — >• {$,,y), with T = ^y, ti = ^(1 — y), and the Jacobian l^^^l = C The area of 
integration changes from r > and ri > to ^ > and < y < 1. After this change of 
variables, the quantities ry and r± will be proportional to ^, rendering the integration over ^ 
trivial. The last integral over y can be calculated in a straightforward way, yielding after some 
manipulations 



X2<0 : M'jllf = MBta]-^ 

Dir]{x2) 

X2> : Mo'j'a? = Mb ia ^- — 
Dir]{x2) 



Dj \ -2/3x2 J 



(91) 



(92) 

where 

r]{x) = (1 + a;/5) A + (1 - x/3) ^2- (93) 
The result for Mq^24* is not the same for X2 < and X2 > 0. This is caused by the propagator 



decomposition (P3|). However, the complete result, with the photon-pole residue included, will 
be independent of the sign of X2- 
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4.1.2 Photon-pole residues 



Now we calculate the photon-pole residues. There are two such contributions, one in each of 
the terms in the decomposition (^). We will indicate these two contributions by M0124 and 
^0124 5 respectively. For Mq124 the contour is closed in the lower half of the complex kQ-pla.ne. 
In that case the photon pole is situated at = \k\ — io, yielding 



0124 



Mf 



a 1 — [3x2 
2^ Di 



2ti t]{x) {1 - V2-nk) 



' f23x + io 

In 

\l-px 



\n{D) + \n{-Di) 



■ (94) 



Here stands for the unit vector in the k direction and VLk indicates the angular variables in 
spherical coordinates (with the polar axis defined along p). For M^ijl the contour is closed 
in the upper half of the complex /cg-plane. The corresponding residue can be obtained from 
Eq. ( P^ by adding an overall minus sign and by substituting [3 —[3 and Di D2 inside the 
square brackets: 



Y,2 



0124 



Mb — 



a 1 — (3x2 



27r Di 



2tt r]{x) {1 - V2-nk) 



In 



-2px + io 
l + px 



-ln(D)+ln(-D2) 



• (95) 



Next one can perform the integration over the azimuthal angle, with the help of the formula 

7^ 1 1 

(96) 



27r 1 - Vi-fik 



This expression is a possible source of collinear divergences, which are regularized by introducing 
the small non-zero fermion masses. In terms of this regularization, |x — Xj| is replaced by 
x — XiY + mf{l — Xi)/Ef. The sum of the two photon-pole residues amounts to 

1 



M0L4 



a 1 — (3x2 



27r Di 



dx 



ri{x) |x — a;2| 



The last integration gives rise to integrals of logarithmic and dilogarithmic type. Let us single 
out the answer for the ^-function-dependent terms: 

a 1 — (3x2 . 



^'^-'0124 



Mb — 



2t\ Dir]{x2] 



in 



{c,[e{-x2) - e{x2)] + 2C2e{x2) + 2C^e{-x2)}, 



(98) 



where 



^ Ki^)-"(^)-<3)' --K^^?iy)-<- 



+ X2 



X2 



7]{- 



X2 



(99) 



Separately, the particle-pole residue and the photon-pole residues depend on the sign of X2. 
However, the sum of these terms does not. This dependence on X2 at the intermediate stage of 
the calculation is a consequence of the decomposition of the unstable VT-boson propagators. 
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The final answer for tlie contribution of the infrared-finite virtual four-point function -D0124 
to the non-factorizable matrix element is given by 



Mr 



0124 



Mr. 



a 1 — 13x2 
2tx Dir]{x2) 



F2{X2\X2) - F2{-D^\X2) 



In 



Do 



2 In 



+ In , \ + 21n, 

\(3{l-X2) 



+ >n(4) 



where 



I D1 + D2 



(100) 



(101) 



(3 D^-D2 

The logarithmic and dilogarithmic functions Fi 2 can be found in App. lCl] . The final answer for 
M0124 agrees with the answer presented in 0. It is also in complete numerical and analytical 
agreement with the corresponding expression in Sect.^ which was calculated with the help of 
the MST. 

The contribution from the other infrared- finite virtual four-point function, -D01235 can be 
obtained from Eq. (|100|) by substituting (pi, ki) ^ {p2, ^2)- 



4.1.3 The pure photon-pole part 



As was already explained in Sect. the photon-pole parts of the virtual scalar functions 
can be related to the corresponding bremsstrahlung interferences. To this end, one needs to 
calculate the pure photon-pole contribution to the matrix element, without performing the 
decomposition of the unstable VT-boson propagators, since this decomposition mixes photon- 
and particle-pole contributions. 

This calculation is pretty much the same as the one discussed in the previous subsection. 
We present only the answer: 



^0^24 



Mb — 



a 1 — (3x2 



2tc Diri{x2) 



m2i 

D2 



F2(X2|X2) - F2i-Do\x2) - fA-D^' I3\x2 



+ Fi [-Do] -f3\x'^ + Fi 1^x2; [3\x-^ - Fi [x2] -^\x^ 



■ (102) 



Note that the photon pole has been evaluated in the upper half of the complex kQ-plane. 
The reason for this lies in the fact that we have opted to perform the calculations in the most 
economic way. In this approach -D0124 is obtained from -D0123 by substituting (pi, ki) (p2! ^2)5 
which automatically shifts the photon-pole from the lower to the upper half-plane (see Sect. |2.4| ). 
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4.2 The infrared-divergent scalar four-point function 



In this subsection we present the calculation of the infrared-divergent virtual scalar four-point 
function -Doi34- In the DMI method such functions are not needed for the calculation of the 
non-factorizable corrections. They arise only in the form of initial-final state interferences. 
Such corrections vanish when the corresponding bremsstrahlung interferences are taken into 
account, as was explained in Sect. |2.4 We perform the calculation mainly to study how one 
can handle infrared and collinear divergences in the DMI scheme and to provide an independent 
check of the results obtained in Sect.^j. 

The infrared-divergent virtual scalar four-point function -D0134 is defined as 

f d^k 1 
" J (2^ [A;2 - A2 + 10] [Di - 2(pi • k)] [-2{ki ■ k) + 10] [2{k2 ■ k) + lo] ' 

We regularize the infrared divergences by introducing a regulator mass A for the photon. There 
are also collinear divergences, which are regularized by the small non-zero fermion masses. 

The pole structure of this integral is such that no propagator decomposition is required. 
There is one photon pole in each of the half-planes of the complex variable k^. There are two 
particle poles in the upper half-plane, and only one in the lower half-plane. Therefore, we opt 
to close the integration contour in the lower half-plane, resulting in only two contributions to 
the scalar function: one photon-pole residue and one particle-pole residue. 



4.2.1 ParticIe-poIe contribution 

One can proceed in the same way as in Sect. [4.1|. We take the residue at the particle pole 
ko = V2-k and exponentiate the propagators by introducing an integration over "time" : 

D^lll = i JdT dt e'-^ I ^ ^ "—^ , (104) 

where f = r {v — V2) + t {vi — V2)- Again we can perform the integration over the momentum 

— * 

k, which is similar to the A-screened relativistic Coulomb potential of a moving particle. As 
the scalar function is infrared- divergent, one should keep the photon mass A. The result of the 
integration is 



lrl+rl{l-vl) 



Here both A and 1 — vf are small. We will consider the limit A — > and Vi ^ 1, such that 
A < Jl-vl 
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The particle-pole contribution now takes the form 



r2+ri(l-^i) 



where + (1 — = a + 6t + ct^, with coefficients 



2 

2/1 ^? \2 I ^2 2 n2 ■ 2q 

r (1 - (3x2) + y ^2, 



777/ 

2r(l - X12) (1 - (3x2) + 2 -| T/3 sin ^2 sin Ou, 

^2 



m 



2 

1 -Xi2)' + ^sin2^i2. (107) 



E2 

The integral over t is logarithmically divergent in A: 



-.\fa+bt+cfi 







-C + ln 



VaTbtTct^ Vc 

where C is the Euler constant. 

The last integration over r is relatively simple, yielding 



XE2ib + 2y/E3] 



T-)part 
-'^0134 



87rsi2 Di 



ln(^Vlnte-lUln^^ 



iMlJ \Ml 'J "'\m2 



(-) 



The invariants S12 and S211' are defined in Eq. (^61). 
4.2.2 Photon-pole contribution 



(108) 



(109) 



Next the photon-pole residue dX = uj = \J k?' + \^ — io is determined: 

d'^ f ^ ^ fix 1 

''''^''16EE,E2J (27r)3u; (1 _ ^a;)[u; - 2^(^][cu - (n,-t;i)][cu - |^| (n^-tJa)] ' 

We want to keep the propagators [u — \k\ (fik-Vi)] as they are, instead of writing them as 
\k\[l — fik-Vi], as was done in Ref. H]. Keeping the exact form of the propagators leads to 
double-logarithmic collinear divergences. If, instead, the simplified version were to be used, 
then the double-logarithmic terms would be lost, and one cannot be sure whether the single- 
logarithmic divergence and the finite part would be correct. 

First we perform the integration over \k\. The presence of A in the light-fermion propagators 
complicates things considerably. The light-fermion propagators can be rewritten in the following 
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way: 



[u - \k\ {nk-vi)]{u} - \k\ {nk-V2)] \k\ [uk-vi - nk-V2) 



uj -\k\{nk-vi) uj -\k\{nk-V2) 



After the integration over \k\ the photon-pole contribution will be of the form 



-'^0134 



WEiE2D] 



[h + h + h 



'112) 



where 



In 



A-n [l-nk-Vi][l-nk-V2] \E\{l-(3x 



h 



V 



Att [rik-vi - nk-V2][l - {uk-Vi)'^] 



In 



1 - Uk-Vi 



;ii3) 



The expression for I2 can be obtained from Ji by the substitution vi ^ V2- The integral Jq 
is similar to the one that shows up in the approach of |Q. It contains only single- logarithmic 
coUinear divergences. The integrals /i^2 are new and give rise to double-logarithmic terms. They 
are evaluated as a principal- value integral, since the singularity present in l/[nk-vi — Uk-vo] is 
an artefact of the split-up ( 111 ) and disappears in the sum Ji + l2- 

As a next step we integrate over the azimuthal angle 0. We obtain for the first integral 



dx 



2(1 - XI2 



Ji — xKi _^ J2 — XK2 



\X — Xi\ 



X — X2\ 



with 



+1 
-/_ 

) {X - Xa){x - Xb) 

Xi + X2 ± i sin 9i sin 62 sin 



In 



EX{l-(3x) 



(114) 



Jl 
J2 



1 + X12 

1 - a;i2 - xi(a;i - 0:2), Ki = X2 - xiXi2-, 

1 - X12 - X2(a;2 - a;i), ^^2 = a^i - 3^23^12- (115) 
As indicated in Sect.[4.1|, |x — Xj| should be regularized by keeping the small non-zero fermion 



masses: \x 



Xi I 



{x — XiY + fif, with /i^ = m^(l — xf)/Ef. Using the result for the 
principal- value integral given in App. |C.2"| , we find for the second integral 

1 



dx 



x+ 



Jb a/_|_ 



ViX 



+ 



1 + ViX 



In 



ViX 



1 + ViX 



;il6) 



where x± = ^^(1 ± Xi2)/2 . 



One is left with a one-dimensional integration over x. The integrals /i^2 can be expressed 
in terms of the dilogarithmic functions J^i and J-" 2, defined in App. |C.3| : 

1 



h,2 



8x_ 



[J^i(x+;fi,2) +J^2{x+)]. 



;il7) 
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The answer for the integral Jq is 



1 ( E\j3\ fmim2\ ^ Xi-Xa ^, 1, 2^ 

where the function /C is introduced in App. |C.3 . 

The complete expression for the infrared- divergent scalar four-point function is obtained as 
the sum of the particle-pole and photon-pole contributions. When comparing all the above 
results with the ones obtained in Sects. |3.2.3| and |3.2.4| , complete numerical agreement is found 



for both virtual and real four-point functions, including coUinear and infrared divergences. 



4.3 Non-factorizable corrections from the five-point functions 



In this subsection we describe the calculation of the scalar five-point functions using the DMI 
method. The contribution of the virtual five-point function to the non-factorizable matrix 
element is given by 



M = zMb 



167ra {ki ■ ^2) 



(27r)4 [P - A2 -f io] [-2{k ■ki) + 10] [2{k ■ ^2) + 10] [Di - 2{k ■ pi)] [D2 + 2{k- P2)] ' 

(119) 



In analogy to Sect.^^, one can perform a decomposition of the unstable l^-boson propagators 



[Di-2{p,-k)][D2 + 2ip2-k)] 



+ 



Di - 2{pi -k) D2 + 2{p2 ■ k) 



D + Ap-k 



(120) 



In this way the matrix element splits into two terms. If the integration contour in the complex 
fco-plane is chosen properly, each term involves one photon-pole contribution and one particle- 
pole contribution. 



4.3.1 Particle-pole residues 



We first calculate the particle-pole residue that contributes to the first term in Eq. ( |120| ). We 



proceed in the usual way, by taking the residue aX kf) = V2-k and subsequently exponentiating 
the propagators. In this case this procedure requires three integrations: 

M2--' = Ms — ^ Idtdrdr, e4&-S-] / ^ ^ .^21) 



I ' J (27r)3 [k-V2Y - - A 

where r2 = tv+ti {v—V2)+t {vi—V2)- The integral over d^k is the same as the one evaluated for 
the calculation of the infrared- divergent four-point function [see Eq. ( p.05| )]. Again, a non-zero 
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photon mass A is needed for the regularization of the infrared divergences. The particle-pole 
residue now amounts to 



^'part- 



where Tn + (1 



-Ml 



ia 1 — Xi2 
~8 



rjf + (1 — Vo) 



a + bt + ct'^, with coefficients 



(122) 



[t(3x2 - ri(l - Px2)Y + -d(^ + ^irf^ sin^ 6; 



El 



ml, 



-2(1 - X12) [Tf3x2 - ri(l - f3x2)] + 2 — f ^(r + n) sin sin 62 

En 



rrio 



:i-xuy + ^2 



sin 6'i2. 



(123) 



Following Sect. [4.2| , first the integration over t is performed, yielding the logarithmically- 
divergent result 



-Mf 



8^2 



-C + ln 



^\fcm2 



\E2 (b + 



(124) 



To linearize b + I^Jac with respect to one of the integration variables, one should make a 
change of variables according to (r, ri) {C,,y), with r = C,y, Ti = ^(1 — y). In this way, the 
integration over ^ can be trivially performed: 



M2P'''*= Me 



la 



dy 



-l + ln 



/ 4m2(l 



3^12 



VAE2(6' + 2vVc') 



+ ln 



Dy + 2D^{l-y) 
AiE 



(125) 

where the coefficients a', b', and c' follow from the coefficients a, 6, and c, by substituting r y 
and ri ^ (1 — y). 

The last integration is technically quite involved, but only gives rise to logarithms. Note 
that one should carefully analyse the infrared and collinear divergences, present in this integral. 
The final answer is formally different for X2 < and X2 > 0- This is the same phenomenon as 
observed in Sect. [1.1| , which can be attributed to the decomposition of the H^-boson propaga- 
tors. The result for the photon-pole residue will compensate this dual behaviour, leading to a 
combined result that is analytically the same for both 0:2 < and X2 > 0. In a similar way 
the second particle-pole residue, corresponding to the second term in Eq. ( |120|) , will formally 
depend on the sign of Xi. 

We will not bore the reader with all the different cases and merely present the answer for 
the case Xi > 0, X2 < 0. Taking into account both particle-pole residues, the final answer reads 



Xi > 0, X2 < : 
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, 2(1-/3x2) , . . 2(1+^ fD\^ 4Px2 , . _ .1 

Dir]{x2) D2r]{x2) \DiJ Dr]{x2) J 

-Mb — Ix2^xi, ms^mi, E2 ^ Ei, ^ D2, P ^ ^pl (126) 
Note that the terms In(^) and In(^) cause the difference with the results presented in Ref. 0. 



j^^'part' ^ _ 




4.3.2 Photon-pole residues 



Next we determine the photon-pole residues. Each of the terms in the propagator decomposition 
gives rise to one photon-pole residue, situated at ko = u = ±\/k'^ + X"^ — io. In the same 



way as in Sect. [4.2| , the light-fermion propagators occurring in the photon-pole residues can be 
rewritten according to Eq. ( p.ll| ). Again we introduce spherical coordinates, with the polar axis 
defined along p. For the integration over we keep the A dependence of u in order to get the 
correct divergences. The combined result of all photon-pole residues is given by 



M' 



with 



- a , 
-Mb - (1 - X12) 

TT 



^{D,,D2,x) 



d^nki ^{D,,D2,x) 



An 



ai)(l-a2) D1D2 



$(«!, 0:2) + $(a2, ai] 



(127) 



D1D2 

1-px 
Dir]{x) 
2(3x 



In 
In 



\E 



+ ITT 



(1 - Px) 



l + [3x 



In 



D2vix) \M^{l + /3x) 



Drj{x) 



iTT [9{Px) - 9{-Px)], 



$(ai,a2) 



{rik-Vi) for {i 



(tti - a2)(l - "i! 
1,2). 



In 



1 — «! 



:i28) 



and ai 

The \&-term in Eq. ( |1271 ) also emerges in the calculations presented in , up to the divergent 
term ln(Ai?). This \E'-term contains infrared divergences and logarithmic collinear divergences. 
The other two terms in Eq. ( |12?| ) are of a type that was already encountered in Sect. ^72[ They 
will give rise to double-logarithmic collinear divergences 



As in Sect. |4.2|, we proceed by performing the azimuthal integration. The remaining inte- 
gration over X = cos 6 gives logarithms and dilogarithms. Most of the ingredients of this final 
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step in the calculation have already been discussed in the previous subsections. Therefore we 
only give the answer. First we do so for the \l/5i-terms. As was observed for the particle-pole 
residues, the results depend on the sign of Xi^2- Adopting the same sign choice as in Eq. ( |126| ), 
we obtain 



xi > 0, X2 < : 



a 



MB — i7T — <n. 



2n D 



V,2 



In 



77(1) 



2 In 



D 



J=a.,l 



In 



1 — Xi 



ln(^U21nr^(^^^ 



J=a.,l 



In 



D 



^ > + 2 In 



1 — Xi 



+ 7^2 

X2 — Xj 
—Xj 

-7^l 

Xi — X 



-Xj 



In 



In 



ml 



The coefficients TZ are given by (i = 1, 2) 

KiD + /3Ji {D, - D2) 



n 



2(3D 



1 + X12) ri{xa) ri{xb) ri{xi)' 



7^i 



2(3x^ 

7]{Xi) 



and IZa t 



+ 2 In 



-X2 
I - X2 



+ 2 In 



Xi 
l+Xi 



I^Xa^i 



(129) 



(130) 



As in Sect. only the sum of the particle-pole residues and Mg^ is independent of the sign 

of Xi o. 



The evaluation of the remaining terms in Eq. ( |127| ) is straightforward. The answers for the 
^^0- and $-terms are given by 



a 1 fmim2\ 
-Mb — „ „ In' ' 



-Mb — 



71 D1D2 \ S12 

a 1 — X12 1 



2tt D1D2 4,T 
The functions J-'i and J-" 2 can be found in App. |C.3| . 

Finally, the answer for the \E'i2 -term reads 



[J^i{x+;v2) +Ti{x+;vi) + 2J^2ix+)]. 



with 



-M, 



B 



ail — I3x2 
2i\D^r^{x2) 



V mi 



In 



-Fi(x2;-/3|x2) 



In 



F2{Xi\x2) - Fi{Xi] -(5\X2) 



R 



77,2 



D 



In 



Di 
Ml 



F2i-Do\x2)-F,i~Do;-P\x2) 



(131) 
(132) 



(133) 



a 



-MB^iD,^D2, p^-py 



(134) 
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Here Dq is defined in Eq. (|101|) and the functions Fi and F2 are given in App. C.l . Note tliat 
tlie coefficient i?^ 2 depends on (3. The contribution can be obtained by substituting 
(£^2, "^2, X2) ^ {El, mi, Xi) in Eq. (pil) . 

The final answer for the contribution of the virtual five-point function to the non-factorizable 
matrix element can be obtained as 



(135) 



with the various contributions given by Eqs. ( |126D and (|129|) -( p!34D . This answer was compared 
numerically with the corresponding MST-expression in Sect.^, which was derived by means 
of a decomposition of the five-point function into a sum of four-point functions. A complete 
numerical agreement was observed. 



4.3.3 Pure photon-pole part 

In order to calculate the real-photon radiative interference corresponding to the five-point func- 
tion, one has to determine the photon-pole residue in the lower half-plane, without performing 
the propagator decomposition. The calculation is more or less the same as the one discussed 
in the previous subsection. 

The answer can be written as 

= M(J + M7 + + M^. (136) 



Note that the contribution is the same as before [see Eq. (132)]. The other contributions 
are changed slightly: 



and 



Rri,2 



-Do\x2)\. (138) 



D 

The contribution can be obtained by substituting {E2, m2, X2) {Ei, mi, xi). 



5 Complete results 

Up to now we have focused on the case of purely leptonic final states. For the purely hadronic 
ones there are many more diagrams, as the photon can interact with all four final-state fermions. 



34 



In order to make efficient use of the results presented in the previous sections, we ffist in- 
troduce some short-hand notations based on the resuhs for the purely leptonic {LL) final 
states. These short-hand notations involve the summation of virtual and real corrections to 
the differential cross-section. For instance, the virtual corrections originating from the ffist 
diagram of Fig. IT] can be combined with the corresponding real-photon correction into the 
contribution dai[{ki;k[\p2). In a similar way, virtual and real five-point corrections can be 
combined into da^^l{ki; k[\k2:, k'2) . The gauge-restoring "Coulomb" contribution will be indi- 
cated by da'^{pi\p2)- In terms of this notation the non-factorizable differential cross-section for 
purely leptonic final states becomes 

daLiiki] k[\k2; k'^) = da^liki, k[\p2) + da'fl{k2; /cabi) + c^o-IlI^i; ^11^2; k2) + da^ {pi\p2) . (139) 

Analogously the non-factorizable differential cross-section for a purely hadronic final state (HH) 
can be written in the following way 



daHH{ki; k[\k2; ^2) = 3 x 3 



1 daiiik^; k[\p2) + I da['lik[; k^ + ^ da^S.{k2^ k'M 



+ f da^^M- fcsbi) + \-\ daf^{k,- k[\k2, k'2) + \-\ dafM; h\k2; k',) 
+ \-\ da'flik,; k\%- k2) + \-\ daf^{k\- k,\k',; fc^) + da''{p^\p2) 



(140) 



In order to keep the notation as uniform as possible, the momenta of the final-state quarks are 
defined along the lines of the purely leptonic case with ki {k'^ corresponding to down- (up-) 
type quarks. If one would like to take into account quark-mixing effects, it suffices to add the 
appropriate squared quark- mixing matrix elements (jV^jp) to the overall factor. Note that top 
quarks do not contribute to the double-pole residues, since the on-shell decay W ^ tb is not 
allowed. Therefore the approximation of massless final-state fermions is still justified. 

For a semileptonic final state (say HL), when the decays hadronically and the W~ 
leptonically, one can write 



daniiki; k[\k2\ k'2) = 3 



i da^^lik,; k[\p2) + I da^^l{k[; kM + ^^^(A;^; k',\p,] 



+ I daf^ikr- k[\k2; k'^) + ^ daf^iK'^ h\k2; k'^) + da''{p^\p2) 



(141) 



Upon integration over the decay angles, the functions d(j^l\ and dafl become symmetric 
under ki ^ k[. As a result, expressions dliOl ) and ([141D take on the form of ([1391) multiplied 
by the colour factors 9 and 3, respectively. These are precisely the colour factors that also 
arise in the Born cross- sect ion. Therefore, after integration over the decay angles, the relative 
non-factorizable correction is the same for all final states. This universality property holds for 
all situations that exhibit the ki k[ symmetry. This means that there is only sensitivity to 
the charges of the decaying particles. 

At this point one may wonder how non-factorizable corrections affect Z-pair-mediated and 
Zif-mediated four-fermion final states. In those cases, only five-point functions contribute. 
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of which there are four contributions, as in Eq. ( |140|) . However, in contrast to Eq. ( |140| ), the 
charge factors are pair-wise opposite, such that integration over the decay angles leads to a 
vanishing result. Thus 0{a) non-factorizable corrections to invariant-mass distributions in 
Z'-pair-mediated or ZH-mediated four-fermion processes vanish. This can be viewed as a 
consequence of the zero charges of the decaying particles. 



6 Conclusions 

In this paper we studied two methods to evaluate non-factorizable QED corrections in the 
double-pole, soft-photon approximation. We derived results for H^-pair production, which are 
valid a few widths above threshold. 

One technique (DMI) is an extension of that of Ref. |^ in the sense that the virtual and 
real photonic corrections are clearly separated and also regularized by a photon mass A and 
charged-fermion masses mi and 1712. The resulting formulae are rather complicated and are 
different from those of [^. 

The second method (MST) extends the standard technique in the sense that five-point 
bremsstrahlung interference terms are decomposed into four-point terms and that the soft- 
photon approximation is used from the start in the evaluation of real and virtual ra-point 
functions. The results obtained with this method are much simpler than the DMI ones, but 
the two are in complete numerical agreement. 

The MST can be easily generalized to more involved final states by a straightforward ex- 
tension of the decomposition of five-point functions to the decomposition of n-point functions. 

The methods and most of the actual formulae in this paper can also be applied to ZZ, ZH 
production and to top-quark pair production with subsequent Wb decays. In the latter case 
the top-quark, W and b take the role of W, v and the gluon that of the photon. In that case 
the DMI formulae are directly applicable since no assumption on the neutrino mass was made. 
The MST formulae would need a small modification. 



A Fey nman- parameter integrals 
A.l The on-shell four-point function 

In this appendix we use Ref. to present a compact expression for the on-shell four-point 



function -D1234, which appears in the decomposition of the virtual five-point function in Sect. 
The result for D^2?,\ can be obtained from -D1234 by the substitutions pi —pi and ki —ki. 
The results of provide an independent check on the formula presented here. 

As was mentioned before, the four-point function D1234 has to be calculated without soft- 
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photon approximation and in the on-shell hmit. The resulting function, defined as 

^ _ [ d^k 1 
'''' ~ I (27r)4 [{k - - Ml.] [{k + p2? - M^] [{k -hY- ml] [{k + k^f - ml] ' 

does not contain any divergences. Upon neglecting the fermion masses mi and m2, the answer 
reads 

M^a(xi - X2)Di2M = ^(-I)''|£i2fri4; -Xk) + -Ciaf— ; -Xk 



16vr^^, [ V J Vri4 



Ml, 'V 2{p^-k 



w 



-^^2 TTo ^o] -Xk - £i2 - —, + io] -Xk + ln(-a;fc) In 



-2 



{h ■ k2) 



(143) 



where ri4 is a solution of the equation 



1 s - 



The quantities Xi^2 are the solutions of the equation 

ax^ + hx + c + ido = 0, (145) 

with coefficients 

a = 2{ki ■ k2) - 2{p2 ■ ki), c = 2{ki ■ k2) - 2{pi ■ k2), 



, ,,2 , 4(pi ■ fc2)(p2 • fei) 4(pi ■P2){ki ■ k2) o(u u \ n/ifi^ 

b = Mw + , d = -2{ki-k2). (146) 



A. 2 The infrared-finite four-point function 

In this appendix we briefiy describe some of the details of the Feynman-parameter integral 
belonging to Dqus (see Sect. ^•2.TD . The integral to be evaluated is given by Eq. (p^. In the 
notation adopted in Sect.j^ this integral reads 







with 



P'(^) ^2^1 _ ^) + -C + Xs + 1/x,) + + 6C + ^(1 - - 6)'. (148) 



AMI, ' ' ' ' 



We perform the following change of variables: 



6 = —^ and ^2 = — (149) 
1 + t + u 1 + t + u 
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Accordingly, the area of integration changes from 

ei>o, e2>o, ei + 6<i, (iso) 

to 

0<t<cx), 0<M<oo. (151) 
The Jacobian of the transformation is given by 



d{tu) 

The final integral to be evaluated looks like 



[1+U + tf 



(152) 



7 7 

M^Doi23 = Ir- I ^ -2 • (153) 

16^^2 I [y, + t\[t^ + t{Xs + l/x,) + 1 + n(l + tC) + ^ 



The second expression in the denominator of the integrand is linear in u up to the small 
term v? m\/ M'^, which is needed to regularize the coUinear divergences of the integral. When 
performing partial fractioning of the integrand, this term should not be neglected; it has to be 
treated small parameter. 

Performing first the integration over u and then the integration over t, we obtain the final 
result (1^) for i5oi23, expressed in terms of logarithms and dilogarithms. 



A. 3 The infrared-divergent four-point function 



Next we present a few steps in the calculation of the Feynman-parameter integral belonging 
to -D0134 (see Sect. |3.2.3| ). The first step involves the integration over momentum space, as 
represented by Eq. (|68D. The contour will be closed in the lower half of the complex fco-plane, 

where only one pole is situated at /cq = 
p{C) direction: 



'0134 



A;^ + — io. We introduce cylindric variables in the 
pdp dz 



Vp' + z' + X' [-1^(01 Vp' + ^' + A^- \ v{0\z + A{i) + io 
The integration over p becomes trivial: 

dz 



(154) 



'0134 



(0 



The last integration simplifies if one uses the representation: 

i d { f dz 



2 • 



(155) 



'0134 



(0 



W\\E{S,)\ d\E{i)\ \ J -\E{S,)\^^+l-\m\z + A{i)/\ + 



10 



(156) 
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Introducing the standard variable transformation t = + 1 ~ z , the final integration can 
be performed, leading to the result given in Eq. (|70D. 



The second stage of the calculation involves the integration over the Feynman parameters. 
We start with Eq. ([701). This time we combine a change of variables, analogous to the one 
utilized in the previous appendix, with a rescaling : 

6 = ^ -7, 6 = -r, n followed by t t, u. (157) 

In this way -D0134 takes the form 

where the definitions of A' and p'^ have now changed. Those quantities are rescaled versions 
of A and . The rescaling is performed in such a way that p'^ changes to 



-rT-pr = t +u H tH u + 1. (159) 

AMy^f mim2 m2 rrii 



In order to linearize this expression with respect to m, one has to introduce one more variable 
transformation 

mi mo 

t = t' + cu, with c=^^, (160) 

which leads to 

, , ^9 =t tu-\ t -\ u + 1. (161) 

4M^ mim2 m2 mi 

After changing the order of integration according to 

0000 0000 0000 Ooo 

J duj dt ^ J du J dt' = J dt' J du+ J dt' J du, (162) 

00 -cu -00 -t'/c 

one can perform the rest of the Feynman-parameter integrations to obtain the final result (|72p 

for L'oi34- 



B Why IZ vanishes 

In this appendix it will be shown that the second term in Eq. (0), given by 

is actually zero. In this integral the photon is not necessarily on-shell, because the residue is not 
taken in the photon pole. However, by power counting we can conclude that only soft photons 
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(ko) 





"^oj* •••• 



(a) 



(b) 



Figure 6: The pole structure of TZ in the soft-photon approximation (a) and the transformation 
of the integration contour in the complex /co-plane (b). The solid circles indicate the particle 
poles, the open circles the photon poles. 

give a noticeable contribution to the integral. All other contributions are formally of higher 
order in the expansion in powers of Tw/M^. Therefore we use the soft-photon approximation 
to evaluate this integral. As a result, all particle poles are situated in the same half-plane of 
the complex ko variable, as is shown in Fig. ^a). 

Next one can deform the integration contour in the way depicted in Fig. ^(b). Note that the 
orientation of the contour is reversed. Figure ^ shows that the sum of the particle-pole residues 
is equal to the sum of the photon-pole residues with the opposite sign. This is a consequence 
of the soft-photon approximation and of the fact that all particle poles turned out to be in the 
same half-plane of the complex variable. The latter is the result of the transformation ( PBD 
introduced in Sect.|2.2|. 

Let us consider the following general integral 



where p^^ is an arbitrary vector. In the soft-photon approximation the denominators can be 
written as Nq = k"^ — + io and Ni = 2{pi ■ k) + pf — mj + io. As mentioned in Sect.p.2|, the 
momenta pi are time-like and have positive energy components, i.e. Ei > \p^\. For simplicity 
we take the photon to be massless, i.e. = 0, but the arguments that follow do not depend 
on this. Deforming the integration contour as described above, and subsequently taking the 




(164) 
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residues in the photon poles, one can write 
d^k {p ■ k) 



{2tt)^2uj NiN2NsN4 



+ 2ni 



u}=—\k\+io 



dPk 



{p-k) 



(27r)4 2tu NiN2N^Ni 



u)=\k\ — io 



In spherical coordinates this takes the form 



(165) 



n{p) 



it: 



+ m 



d\k\<pnk \kWE + f{nk)\p\) 



1=1 



\kWE-f{n,)\p\) 



d\k\ d^Qk 

(2^)' n \2\k\ {E, - mk)m +Pl-mj + 



1=1 



(166) 



In the second term one can make a change of variables according to \k\ 
to obtain 



■\k\ and 



-nk, 



TZ{p) = in J 



ikmE+fiQkM) 



d\k\ d^^k 

(2'^)" n \-2\k\ (E, + Mn^mi) +p}- m2 + 10 



1=1 



(167) 



This integral is ultraviolet-finite and all poles are situated in the same half-plane of the complex 
variable since Ei > \p^\ and \fi{Qk)\ < 1- By closing the contour in the opposite half-plane, 
one finds TZ{p) = 0. From this it trivially follows that TZ = 0. 



C Special functions and integrals in the DMI method 

C.l The functions Fi and F2 

In this appendix we present the functions Ei and F2, which are used in the calculations in 
Sect. HI. The function Fi is defined as 



Fi{a;p\xi) 



dx 



X — a 



ln(l + px) 



9{x — Xi) — 9{xi — x) 



(168) 



Here a is a complex number with a non-zero imaginary part, and f3 and Xi are real numbers 
with absolute value smaller than 1. The analytical expression for this function is given by 



Fi{a;(3\xi) = -2Li2(-r- — -]+Li2 



+ In 



1 + ap 

(3{a-l) 
l + a(3 . 



1 + a/3 



Li2 



1 + ap 



\ 1 + ap J 



Hl + (3)-2\n{^^^) ln(l + x./3). 



(169) 
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In addition we need this function in the special case a = Xi, without non-zero imaginary part. 
There, the integral is logarithmically divergent. This is a collinear divergence and 

should be regularized by keeping the small non-zero fermion masses. The answer in this case is 

F,(.,n.,) = Hl^^.mn{^)-u{3l±^)-u{0^). (170) 



The other function, F2, is defined as 





f dx 


\Xi) = ^ 






' X — a 



[x — Xi) — 0{xi — x) 



(171) 



For a and x-i the same restrictions as indicated for the function Fi apply. The corresponding 
analytical expressions are 

F2{a\xi) = -2 ln{xi - a) + ln(-l - a) + ln(l - a), (172) 

and 



F2(x,|x,) =ln(^^). (173) 



C.2 The azimuthal principal- value integral 

In this appendix we present the result for the azimuthal principal- value integral, used in Sect.^: 



2n 



'0 



with 

y4 = (fi — t>2 cos6'i2) cos6' and i? = t>2 sin 6*12 sin 6*. (175) 
The principal-value integration yields 



1 



1 



for A/B G (-1-1, +CX)) or equivalently cos6' G (-|-a;_|_, -|-1) 
for A/B E {—00, —1) or equivalently cos6' G (—1, —x^) 



(176) 



where 



X-\ 



\v1+V^- 2Wif2Xi2 



and - 52 = - Jx^ - x\. (177) 



In non-coUinear situations one can take v\^2 1, resulting in x+ = y^(l -(- xi2)/2 . 



42 



C.3 The functions T2 and /C 



In this appendix we present the functions JFi, JF2, and /C, used in Sect.^ for the infrared- 
divergent four- and five-point functions. 



The function T\ is defined as 

1 



dx 



1 /I — vx 
In 



2 1 — vx \1 + vx 



(178) 



Here, x+ is real with < x+ < 1, and the quantity v is real and close to unity. For v ^ 1 the 
answer for this integral is given by 



J^i{x+;v) 



1 



1-x 



2 V 2 



!^ + 1 ln2(l _ xl) - ln(x+) ln(l - ^ 



;i79) 



The function is defined as 
which amounts to 

J^oix 



dx 



x+ 



1 , /l-x 
In 



x^ — X 



2 1 + X Vl + X/' 



2yx+) 



1-x 



LUxt 



"6 



+ ln(x+) ln(l — x\ 



(180) 



(181) 



In our explicit formulae, the functions T\ and T2 always enter as a sum. This sum can be 
represented in a compact form: 



Tx{Xj^\ v) + J'2(Xh 



1 



1-x 



The function /C is defined as 

/C(A;5|xo;/x') = 1 



-lln^fl^ 
2 V 2 



dx 



Li2 



TT 



(182) 



1 (x + A)y^(a;-xo)2 + /i2 



ln(S - x) 



;i83) 



A being a complex number with a non-zero imaginary part, and B being real and larger than 1. 
The quantities Xq and /i are real, with |xo| < 1 and /i^ <^ 1. The resulting analytical expression 
is somewhat more complicated: 



1 f / A+l\ / A + xo\ (B-x^ 
£i2 1;— -£12 1;— -£12 



5 + 1 A + xo 



--In^ 



B — xq 

B-xq A- \ ) 2"' \ B-1 
2 A-1 ' 



B-1' A + xo 
+ \n{B - 1) In 



+ \n{B - xo) In 



4(l-xg) A + xo 



A + 1 
A + xo 

(184) 
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